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We describe an efficient theoretical criterion, suitable for indistinguishable parti- 
cles to quantify the quantum correlations of any pure two-fermion state, based on 
the Slater rank concept. It represents the natural generalization of the linear entropy 
used to treat quantum entanglement in systems of non-identical particles. Such a 
criterion is here applied to an electron-electron scattering in a two-dimensional sys- 
tem in order to perform a quantitative evaluation of the entanglement dynamics for 
various spin configurations and to compare the linear entropy with alternative ap- 
proaches. Our numerical results show the dependence of the entanglement evolution 
upon the initial state of the system and its spin components. The differences with 
previous analyses accomplished by using the von Neumann entropy are discussed. 
The evaluation of the entanglement dynamics in terms of the linear entropy results 
to be much less demanding from the computational point of view, not requiring the 
diagonalization of the density matrix. 
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I. INTRODUCTION 

The entanglement, possibly the most remarkable feature of quantum mechanics, repre- 
sents a fundamental resource for quantum information processing , and the con- 
cept of bipartite and multipartite entanglement is nowadays well-stated for quantum sys- 
tems composed of distinguishable constituents. For an entangled system it is impossible 
to factor its state in a product of independent states describing its parts. On the other 
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land the notion of entanglement is more controversial in systems of identical particles 
^, [5, [g, 0, [s, y, 0, 11 1 . The entanglement of such systems is investigated at present in 



many areas of physics, like quantum 



ultracold gases 



m Q, Q, 



16 



13, 



optics, quantum charge transport in semiconductor, 



18| . Here the difficulties appear in the definition of 



a criterion able to classify and quantify the entanglement. They are mainly due to the ex- 
change symmetry which requires the antisymmetrization or symmetrization of the quantum 
wavefunctions describing fermions or bosons, respectively. 

Different methods to treat the quantum entanglement in systems of indistinguishable 
particles are present in literature. In the approach developed by Wiseman and Vaccaro [l^ 
the entanglement of the particles is a sort of accessible entanglement, i.e. the maximum value 
of the entanglement that could be extracted from the system and placed in quantum registers, 
from which it could be used to perform quantum information processing. In the theory 
introduced by Zanardi the entanglement should be evaluated by using the density matrix in 
a mode-occupation represention and is based on the formal mapping of the Fock space into 
states of qubits [3]. The method proposed by Schliemann is based on the Slater rank of the 
state (i.e. the minimum number of Slater determinants needed to express it) as a counterpart 
of the Schmidt rank criterion usually adopted for distinguishable particles [191, |20| . This 
latest approach has been recently reexamined by other authors, which, for the case of a 
two-particle pure state, have suggested to evaluate quantitatively the entan glem ent as the 
von Neumann entropy (vNE) of the one- p„«e .educed den.t. .a.. BH" 



12J. Here 

the quantum correlations due to symmetrization or antisymmetrization of the wavefunction 
do not represent a genuine manifestation of quantum entanglement j^]. Therefore from 
this point of view such a criterion results to be the natural generalization of the approach 
commonly used to treat quantum correlations in systems of distinguishable particles. 

Following the basic concepts of the approach proposed by Schliemann [J] in this paper we 
discuss an entanglement criterion for two-fermion systems. It is still based on the analogous 
of the Schmidt decomposition theorem for the pure fermion state but requires the calculation 
of the linear entropy (LE) of the one-particle reduced density matrix. We intend to establish 
whether the LE can be considered a valid measure of the lack of knowledge about the 
quantum state describing the system, which should include not only the uncertainty due 
to the impossibility of attributing a definite state to each particle but also the amount of 
uncertainty deriving from the indistinguishability of the particles. In order to get a better 
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understanding of this criterion, we first apply it to study a simple prototype theoretical 
model, then we analyze a system of physical interest, namely a two-electron scattering event 
in a 2D semiconductor structure. The ent ang lement dynamics in such a system has been 
recently investigated in terms of the vNE 15|, l2J] and a comparison with previous results 
lead to the conclusion that the LE can be an efficient and still valid entanglement measure 
for binary collisions, as for other physical phenomena of interest in quantum-information 
processing where identical particles are involved. 

The paper is organized as follows. In Sec. [Ill we describe the main properties of the LE 
as an entanglement measure for two-fermion systems. In Sec. IIIII we evaluate numerically 
the time evolution of the LE for a scattering event between a free propagating and a bound 
electron in a two-dimensional system considering different spin configurations. Conclusions 
are drawn in Sec. [IVl 



II. THE THEORETICAL CRITERION AND ITS EVALUATION IN A 

2iV-MODES SYSTEM 

In this section we introduce the entanglement criterion for two-fermion systems based on 
the concept of LE, usually applied to distinguishable particles. Furthermore we compare 
quantitatively such a criterion with the one based on the vNE in the case of a simple system 
with 2N degrees of freedom. 



A pure state of two fermions can be written as 



19|: 



2N 



\^p)=Y,uj,,a\a]\Q) (1) 

where and a\ are the annihilation and creation operators of the mode i satisfying the 
usual fermionic (anti) commutation rules {aj,aj} = and |0) is the vacuum state. 
are the elements of a complex and antisymmetric {2N x 2A^) matrix VL where 2N is the 
total number of modes for each single particle, while the normalization condition is given 
by Tr[fi^i7] = L The single-particle reduced density matrix p for the state \^ f) can be 
computed from the two-particle density matrix pp = \'^p){'^f\ and its elements are [Qj: 

= J'^'fy^ , = (2) 
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The eigenvalues of p are while the coefficients Zi stem from the Schmidt decomposition 
of \^ f) in terms of Slater determinants Q]. Furthermore it should be noticed that the 
eigenvalues of the one-particle reduced density matrix are pairwise identical and therefore 
it holds \z2k\^ = |-22fc-iP with 1 < A; < A^. The number of coefficients Zk that are different 
from zero is the so-called Slater rank, which can be related to the entanglement as follows: 
a state with Slater rank equal to 1 (i.e. that can be written as a single Slater determinant) 
is non-entangled, a state with Slater rank greater than 1 is a linear combination of two or 
more Slater determinants, therefore it can be considered entangled. 

Many alternative ways of defining a function apt at evaluating the lack of knowledge 
about a subsystem have been proposed in the literature among them the Tsallis entropy 



25 



as 



generalizes the concept of the vNE, encompassing, among the others LE. It is defined 



m ■ 

= ^Tr{p - pn (3) 
q-1 

where g is a real, not necessary positive, number. In the case of q tending to 1 one obtains the 
well known vNE e„jv = — Trjplnp} which satisfies some standard properties as concativity, 
additivity and sub-additivity, and which is acknowledged to be a good quantum correlation 
measure of a pure two-fermion state. When q is equal to 2 Eq. ([3]) reduces to = 1 — 
Trp^, that is the LE 27|]. In this paper we focus on such a quantity as a measure of the 
entanglement in systems of indistinguishable particles. Even if LE is not additive in the 
usual sense as shown in Ref. j23] , it has some interesting properties so far not fully exploited. 
In fact it turns out to be extremely valuable for the application of numerical methods, as 
detailed in the following. 

In terms of LE the quantum entanglement of the pure two-fermion state {"^f) defined in 
Eq. ([1]) is given by 



2N 



2N 



(4) 



From the above expression we observe that the evaluation of the LE can be performed 
directly from the matrix Q thus , in numerical calculations, the definition and allocation of 
the one-particle density matrix p are not required: therefore e^, is much easier to calculate 
than the vNE since no diagonalization of the matrix p is needed. This aspect appears to be 
relevant since the complexity of many systems of physical interest practically prevents the 
diagonalization of the corresponding density matrix. Furthermore we note that Eq. (jl]) is 
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represention-independent since the trace is invariant with respect to unitary single-particle 
t r ansf ormat ions . 

By using the above trace operations we can also express in terms of the eigenvalues 
Izil"^ of the one-particle density matrix p in order to compare the expression of the LE with 
the one of the vNE. By taking into account the above mentioned property for l^;^^, we find 
for the former that 

N 

e, = l-J2'^\zk\' (5) 
while, as shown in the literature j^, [21I, 22], the latter can be written as 



TV 



£„^ = ln2-^2|zfc|2ln2|2;fc|2. (6) 

k 

We stress that Eq. ([5]) has been reported in oder to make more explicit the following dis- 
cussion, but it is not employed in the numerical calculations where Eq. is used instead, 
not requiring the calculations of the Q eigenvalues. In spite of similarities between the two 
expressions for a two-fermion system, and £„jv have a different dependence upon the co- 
efficients Zk- Both Eqs. ([5]) and ([6]) attain their minimum value when the state I^&f) can be 
written in terms of a single Slater determinant. In this case we have {zi]"^ = | while all the 
other coefficients are zero and therefore = | and e^N = In 2. The fact that the minimum 
value of the two measures is not 0, differently from what happens for the distinguishable par- 
ticles case, is related to the unavoidable correlations due to the exchange symmetry. Since 
the quantum correlations related only to antisymmetrization of the state of two fermions 
cannot be used to violate Bell's inequality and are not a resource for quantum-information 



processing (as shown in previous works 2l|, |23|) a state with Slater rank equal to 1 can be 
considered as non-entangled . As a consequence we will assume that a value = | indicates 
a non-entangled state. 

For a maximum correlated state it holds \zk\'^ = W k and = 1 — j^. We note that 
in the case of a two-fermion system with a very large number 2N of modes, the maximum 
value of the LE tends asymptotically to 1 with a power law as for distinguishable particles. 

In order to compare some properties of the LE and vNE here we shall analyze the two 
entanglement criteria for a simple 2A^-modes two fermion system in a state \x), that can be 
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expressed by Eq. ([T]), with the following coefficients of the antisymmetric matrix Q: 



l+{N-l){l-a) 
2N 



2 



for i = 1, j = 2 



for i = 2k- 1, J = 2k with 2<k<N (7) 



2N 

otherwise 



where a is a real parameter ranging between and 1. In particular we observe that for a = 
the elements of the matrix vanish except for which reduces to In this case the state 
\x) can be set in terms of a single Slater determinant and therefore is non-entangled. On 
the other hand when a = 1 the condition of maximum entanglement is reached, with all the 
non vanishing coefficients equal to -^Z^- 

In the above model the one-particle reduced density matrix p is simply a diagonal matrix 
with eigenvalues tu^g c^l^.i 2k- This makes straightforward the calculation of the matrix 
p^. Its trace, needed for the evaluation of the LE , can be written as function of the parameter 
a as: 

2N 2N ^ ^ 

As expected for a = 0, takes its minimum value |, while the maximum value of 1 — 2^ 

is reached for a = 1. The expression of £^ normalized to 1 reads 

el = l-{l- a)\ (9) 
The vNE can also be easily calculated and its normalized form is 



A^lnA^V N 
+ (l + (l-.m-l))ln ^^'^-f^-^) ). (10) 

We note that the normalization of the vNE and of the LE allows us to compare quantitatively 
the two measures of the quantum entanglement. In Fig. [1] we report the entanglement of 
the system as a function of the real parameter a. Both curves get their minimum value 
for a = 0, both increase with a and reach 1 when a gets to 1. From the comparison we note 
that the LE is always greater than the vNE as a consequence of the different normalization 
procedures apart from the initial and final values when they coincide. For a tending to 1 
(i.e. for \x) tending to the maximally entangled state) the two measures attain the same 
value as for the case of distinguishable particles |28|. For the sake of completeness we show 



the two non-normalized curves in the inset of Fig. [H 
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FIG. 1: (Color online) Comparison between the normalized values of the vNE and the LE as a 
function of the real parameter a in the 2N modes two-fermion system described in the text (for 
this calculation N has been taken equal to 2601). The inset displays the dependence of the two 
not normalized measures upon a. 

III. LINEAR ENTROPY DYNAMICS IN A TWO-ELECTRON SCATTERING 

A. The model 



The dynamics of quantum entanglement has been investigated in various physic al p he 



29|, binary collision events 



15 



3 



3l| 



nomena including, for example, ionization processes 

and phonon-atom interaction [s^]. Due to the increasing quest for quantum computing 
capable solid state devices, the study of the entanglement formation in a scattering event 
in semiconductor structures plays an important role. A numerical analysis of the entangle- 
ment dynamics in terms of vNE for a two-electron collision in two-dimensional semiconductor 
nanostructures (GaAs) has been recently presented in literature 15|, |2J] analyzing the evo- 
lution of the entanglement when an electron freely propagating interacts through a screened 
Coulomb potential with another electron bound to a specific site by a harmonic potential. 
In this section we intend to study the entanglement for such a model by using the criterion 
based on the LE described in the previous section. 

The Hamiltonian describing the physical system can be written as: 



+ —] + 

2m V drl drl 



+ -mu (r„ 



rof + -muj^n - (11) 
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where e and m are the GaAs dielectric constant and effective mass, respectively, and tq is 
the center of the harmonic potential, with energy-level spacing hu. Spin-orbit effects have 
not been considered. At the initial time to one of the two particles, namely the incoming 
electron, is represented by a minimum uncertainty wave-packet centered in ro(see Fig. [Jti): 



where a is the mean spatial dispersion, k = y/2rnE^/h with m the effective mass of the 
carrier and Ek is the carrier kinetic energy. The bound electron is in the ground state of a 
two-dimensional harmonic oscillator 



where ri is the center of the harmonic potential, with energy-level spacing hu. The distance 
|ri — ro| is such that at the initial time to the Coulomb energy is negligible. We stress 
that, as a consequence of the scattering, the state of the particle confined in the harmonic 
oscillator (HO) changes and that the two-particle wavefunction can be expressed as a single 
Slater determinant \l/(r, r') = iplr) (f){r') — (f){r) ilj{r'), only at the initial time. In order to 
better analyze the state of the bound particle. Figs. [2)D-d display the square modulus of the 
projection of the antisymmetrized two-particle wavefunction \I^(r, r') at time t = 480 fs on the 
first three eigenstates ^nin = 1, 2, 3) of the harmonic oscillator: 7ra(r) = | J (ir',^„(r')\I^(r, r')p. 
Note that the spectral decomposition on the HO ground state at time t = (Fig. [2^), 
representing the square modulus of the one-particle initial wavefunction of the free carrier, 
has been almost totally transmitted with no reflected part at t = 480 fs . The scattering 
event leaves the HO in a superposition of excited states, as can been seen in Figs. ^ and 

and for higher energies the peaks of the function 7„(r) are closer to the center of the 
HO. This is due to the different energies of the outgoing particle as the bound particle is 
left in different HO excited states. Note that the Coulomb potential also creates spatial 
correlations between the two electrons. 

We consider now the effect of different initial spin configurations on the evolution of the 
entanglement. In the first quantum state studied the two electrons have the same spin (spin 
up): 




(12) 




(13) 




(14) 
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FIG. 2: (Color online) Square modulus of the projection of the antisymmetrized two-particle 
wavefunction on the first three energy eigenstates of the harmonic oscillator centered in ro=(95 
nm, 95 nm): 7n(r) = |/ dr'^„(r')\E'(r, r')p. The two upper graphs show the projection on the 
ground state at two different times t = (a) and t = 480 fs (b), while in the lower graphs the 
projections on first (c) and the second (d) excited state at time t = 480 fs are reported. 

where the wavefunctions corresponding to the states {ip) and |0) are of the type defined in 
Eqs. ( IT2l) and ( fT3l) . respectively, and the ket |t) indicates a spin up state. The explicit form 
of the matrix Q for the state can be obtained by discretizinff the spatial coordinate r 

n 

into a A^^ points grid (A^ points for each spatial dimension) as shown in Ref. [15|]. We get 



for the 2N'^ x 2N'^ matrix 



1 / , 

^^ = ^{ ' (15) 

V2 \^ 

where Qa is the antisymmetric N"^ x A^^ matrix whose elements read 

ujij = ipivi) (j){rj) - 0(ri) ipirj). (16) 

The second case considered is the one with two electrons having different spins, that 
cannot be factorized in a spin and a real space term. The form of the two-particle state is 

|T) = ^(i^0)in)-i</.^)iiT)). (17) 
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By applying a procedure analogous 
transformation for the spin variables 



;o the one used for |\&) and introducing an unitary 
isl . we get 



(18) 




where Qs is the symmetric counterpart of Qa- 

The last two states considered still describe two electrons with different spins, but contrary 
to the previous case they can be factorized in a position term and in a spin term. We can 
identify the singlet spin state 



1$) 



1 



\^4>) + \4>^)][\n)-m 



with 




(19) 



(20) 



and the triplet spin state 



in) + iiT) 



with 




(21) 



(22) 



The L 



15 



of the triplet state can be easily obtained from the one of the same-spin state |\E') 
2l|. In fact for the eigenvalues of the one-particle reduced density matrix of the state 
) it holds 



1 



for l<i< N^. 



(23) 



Therefore from the Eq. ([5]) its LE is given by 

Ar2 



2(1 + - 



(24) 



The property expressed by the relation is due to the fact that the triplet state, like the 
singlet state, cannot be factorized in a space term and a spin term and this property remains 
true during the time evolution. From Eq. 0241) we observe that the minimum value of the 
LE of I S) is I , greater than the value of | corresponding to a non-entanglement condition 
for any two-fermion state according the criterion introduced in the previous section. Such 
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an offset is related to the lack of knowledge about the spin of the particles. The behavior 
shown by the triplet and singlet states is in agreement with the one obtained in previous 
works where the entanglement formation is evaluated in terms of the vNE for the case of 



one-dimensional and two-dimensional scattering 



24| . as we show in the following. 



B. Numerical results 

In order to calculate the entanglement dynamics in the system described above, we solve 
numerically the time-dependent Schrondiger equation for the two-particle wavefunction con- 
sidering as initial condition two electrons described by the wave-packets given in the Eqs. f|T2|) 
and (fT3|l . In this way at each time step we have the two-particle wavefunction needed to 
define the matrix Q. Finally from expression (jl]) and by using only the matrix elements Uij 
we can compute the time evolution of the entanglement in terms of the LE . We stress that 
such an approach is simple and does not require matrix diagonalization procedures which 
result to be very demanding from the point of view of the numerical calculation. 

Figure [3] shows that at initial time the LE for the states and |T) get its minimum 
value |. In fact the Coulomb energy is still negligible being the two wave-packets far enough, 
and the only quantum correlations present are due to the exchange symmetry. Therefore, 
as expected, and |T) must be considered as initially non entangled. We observe that as 
the free carrier get closer to the center of the harmonic potential the quantum correlation 
builds up and the LE reaches a stationary value once the particles get far enough. Such a 
value depends upon the initial kinetic energy of the propagating carrier Ek: in particular 
it is higher for higher energies for both spin configurations, in good qualitative agreement 
with the previous results found for scattering events between two distinguishable particles 



33, 



34 



35| 



In Fig. mthe time evolution of the entanglement for the singlet |$) and triplet |S) spin 
states are presented. As expected, at the initial time the LE is |. This implies that |$) and 
|S) are initially entangled, being their Slater rank greater than 1. In fact from Eqs. (fT9|) 
and fl2Tl) we observe that they cannot be put in terms of a single Slater determinant. Such 
a result confirms the correspondence between the LE and the Slater rank criterion. 

In order to better compare the properties of the entanglement measure for an electron- 
electron scattering obtained using the LE and vNE, the time evolution of the entanglement, 
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Same-spin Ej,= 10 
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FIG. 3: (Color online) Entanglement as a function of the time for different initial states: same-spin 
state 1^) for two different initial energy values of the incoming electron = 10 meV (dashed 
line) and 20 meV (dotted line); the state with electron having different spin |T) for = 10 meV 
(dot-dashed line) and 20 meV (solid line). The harmonic oscillator energy is fko = 2 meV. 
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FIG. 4: (Color online) Entanglement as a function of the time: triplet spin state |H) for two 
different initial energy values of the incoming electron Ei^. = 10 meV (dashed line) and 20 meV 
(dotted line) and singlet spin state 1$) for E^ = 10 meV (dot-dashed line) and 20 meV (solid line). 
The harmonic oscillator energy is /iw = 2 meV. 
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FIG. 5: (Color online) Two measures of the entanglement as a function of the time. The time 
evolution of the normalized LE for same spin state |^') (dashed line) and for the state with electrons 
having different spin |T) (dotted line), is compared with the time evolution of the normalized vNE 
for the same states: l^') (solid line) and |T) (dash-dotted line). The initial kinetic energy of the 
incoming electron is 30 meV while the harmonic oscillator energy is ^ = 2 meV. 

evaluated according to the two criteria, is presented in Fig. [5] for the states and |T). As 
for the case of the theoretical model studied in the previous section, here we have normalized 
the two measures in order to compare them quantitatively. At initial time both of them 
are zero since no quantum correlation is initially present apart from the one related to the 
exchange symmetry. Then, at increasing times, for a given state the LE is always greater 
than the vNE. Nevertheless the time of the entanglement formation, defined as the time at 
which the entanglement reaches its stationary value, is the same for both measures. This 
behavior can be ascribed to the fact that, as indicated by the Eqs. ([5]) and ([6]), both the LE 
and vNE can be expressed as function of \zk\'^, the eigenvalues of the one-particle reduced 
density matrix, which can be assumed weakly time-dependent for large times. 



In the last years the notion of the entanglement for systems of identical particles has been 
widely discussed: various approaches have been proposed in literature each having advan- 
tages and drawbacks 0, Q, Q. In this paper we have analyzed one possible theoretical 



IV. CONCLUSION 
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criterion to quantify the quantum correlations appearing in a two-fermion systems. It can 
be considered the generahzation of the LE, usually adopted in the context of distinguishable 
particles. Such a criterion is based on the fermionic analogous of the Schmidt decomposition 
theorem and uses the one-particle reduced density matrix. In particular our analysis shows 
that LE permits to determine whether the uncertainty concerning the states derives only 
from the indistinguishability of the particles or it is a genuine manifestation of the entan- 
glement. This aspect is crucial since, as shown in previous works, the quantum correlation 
due to the exchange symmetry does not represent a good resource for quantum information 



processmg 



21 



23|. 



The criterion proposed here appears to be closely related to the one involving the evalu 



22| . To compare the two criteria 



ation of the vNE of the reduced statistical operator 
and analyze their properties we have quantified the entanglement according to the two mea- 
sures in a simple theoretical model describing a two-fermion system with 2N degrees of 
freedom. We found that the LE is greater than the vNE, but the two measures give the 
same value when the quantum state of the system tends to the maximally entangled one in 
agreement with results obtained for distinguishable particles 281 ]. 

Furthermore we have used the LE to quantify the time evolution of quantum correlations 
in a numerically simulated electron-electron collision event. In agreement with the previous 



analyses obtained by using the vNE 



16|, our numerical results show that the entan- 



glement dynamics depends on the spin components of the states even if the Hamiltonian 
does not include spin-terms. Moreover also in this case the triplet and singlet spin states 
are initially entangled and such an entanglement can be ascribed to the lack on knowledge 
about the spin state of a particle in a specific real-space state. For our two-electron model 
the values of the entanglement obtained by using the LE results to be higher than the ones 
obtained by using the vNE in agreement with what found analytically for a simple model of 
a 2A^-mode two-fermion system. Most notably we found that the time of the entanglement 
formation is the same for the two measures. 

Finally we note that the calculation of the entanglement in terms of the LE is easier 
and computationally much faster than the one performed by means of the vNE since no 
diagonalization of the one-particle reduced density matrix is required. Therefore we believe 
that the LE is an useful correlation measure for a two-fermion system and its application 
turns out to be very helpful to investigate the entanglement dynamics for those physical 
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systems with a very large number of degrees of freedom, whose complexity does not allow 
the diagonalization of the reduced statistical operators through numerical procedures. 
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